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We present a measurement of the forward-backward charge asymmetry (Apb) in pp -^ Z/j* + 
X -^ e+e^ + X events at a center-of-mass energy of 1.96 TeV using 1.1 fb^^ of data collected 
with the DO detector at the Fermilab Tevatron collider. Afb is measured as a function of the 
invariant mass of the electron-positron pair, and found to be consistent with the standard model 
prediction. We use the Afb measurement to extract the effective weak mixing angle sin^ 9'^ = 
0.2326 ± 0.0018 (stat.) ± 0.0006 (syst.). 

PACS numbers: 13.85.-t, 13.38. Dg, 12.15.Mm, 12.38.Qk 



In the standard model (SM), the neutral-current cou- 
plings of the Z bosons to fermions (/) at tree level are 
defined as 



2 cos dw 
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where 9w is the weak mixing angle, and gy and gy^ 
are the vector and axial-vector couplings with gy = 
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2Qf sin 6\Y and g^ = I. 
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Here /g is the weak 
isospin component of the fermion and Qf its charge. 
The presence of both vector and axial-vector couplings 
in qq —> Z/j* —^ i~^£~ gives rise to an asymmetry in the 
polar angle (9) of the negatively charged lepton momen- 
tum relative to the incoming quark momentum in the 
rest frame of the lepton pair. The angular differential 
cross section can be written as 



da 
dcos9 



A{l + cos'^9)+Bcoa9, 



(2) 



where A and B are functions dependent on /g , Qf, and 
sin^ 9w Events with cos6' > are called forward events, 
and those with cos 61 < are called backward events. 

The forward-backward charge asymmetry, Aps , is de- 
fined as 



A 



FB 



ap — <7b 
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where apm is the integral cross section in the for- 
ward/backward configuration. We measure Aps as a 
function of the invariant mass of the lepton pair. To 
minimize the effect of the unknown transverse momenta 
of the incoming quarks in the measurement of the for- 
ward and backward cross sections, we use 9 calculated in 
the CoUins-Soper reference frame [l]. In this frame, the 
polar axis is defined as the bisector of the proton beam 
momentum and the negative of the anti- proton beam mo- 
mentum when they are boosted into the rest frame of the 
lepton pair. 

The forward-backward asymmetry is sensitive to 
sin^ 0^, which is an effective parameter that includes 
higher order corrections. The current world average value 
of sin^ 9°^ at the Z-pole is 0.23149 ± 0.00013 gj- Two 
sin^ 9^ measurements are more than two standard de- 



that from the 

0,f) 



viations from the world average value 
charge asymmetry for b quark production {A^p°g) from 
the LEP and SLD collaborations [3| and that from neu- 
trino and antineutrino cross sections from the NuTeV 
collaboration [^. The Ap^ measurement is sensitive to 
the couplings of h quarks to the Z boson, and the NuTeV 
measurement is sensitive to the couplings of u and d 
quarks to the Z boson, as is the measurement presented 
here. Previous direct measurements of u and d quark 
couphngs to the Z are of limited precision [a, |6|. Thus, 
modifications to the SM that would affect only u and d 
couplings are poorly constrained. In addition, ApB mea- 
surements at the Tevatron can be performed up to values 



of the dilepton mass exceeding those achieved at LEP 
and SLC, therefore becoming sensitive to possible new 
physics effects [7|, |8| . Although direct searches for these 
new phenomena in the Z/^* -^ l^(.~ final state have 
been recently performed by the CDF and DO collabora- 
tions [9|, charge asymmetry measurements are sensitive 
to different combination of couplings, and can provide 
complementary information [lO|. 

The CDF collaboration measured Aps using 108 pb~^ 
of data in Run I [ll| and 72 pb^^ of data in Run II [^. 
This analysis uses 1066 ib 65 pb~^ of data I4I collected 



with the DO detector [13| at the Fermilab Tevatron col- 
lider at a center-of-mass energy of 1.96 TeV to measure 
the ApB distribution and extract sin^ 9'^. 

To select Z/"f* events, we require two isolated elec- 
tromagnetic (EM) clusters that have shower shapes con- 
sistent with that of an electron. EM candidates are 
required to have transverse momentum px > 25 GeV. 
The dielectron pair must have a reconstructed invariant 
mass 50 < Mee < 500 GeV. If an event has both its 
EM candidates in the central calorimeter (CC events), 
each must be spatially matched to a reconstructed track 
in the tracking system. Because the tracking efficiency 
decreases with magnitude of the rapidity in the end 
calorimeter, events with one candidate in the central and 
one candidate in the end calorimeter (CE events) are re- 
quired to have a matching track only for that in the cen- 
tral calorimeter. For CC events, the two candidates are 
further required to have opposite charges. For CE events, 
the determination of forward or backward is made ac- 
cording to the charge of the EM candidate in the central 
calorimeter. A total of 35,626 events remain after appli- 
cation of all selection criteria, with 16,736 CC events and 
18,890 CE events. The selection efficiencies are measured 
using Z/j* -^ ee data with the tag- probe method 14| . 
and no differences between forward and backward events 
are observed. 

The asymmetry is measured in 14 Mge bins within the 
50 < Mee < 500 GeV range. The bin widths are de- 
termined by the mass resolution, of order (3 — 4)%, and 
event statistics. 

Monte Carlo (MC) samples for the Z/7* -> e+e" 
process are generated using the pythia event genera- 
tor [15] using the CTEQ6L1 parton distribution functions 
(PDFs) [16], followed by a detailed GEANT-based simula- 
tion of the DO detector 17]. To improve the agreement 
between data and simulation, selection efficiencies deter- 
mined by the MC are corrected to corresponding values 
measured in the data. Furthermore, the simulation is 
tuned to reproduce the calorimeter energy scale and res- 
olution, as well as the distributions of the instanteneous 
luminosity and z position of the event primary vertex 
observed in data. Next-to-leading order (NLO) quantum 
chromodynamics (QCD) corrections for Z/'y* boson pro- 
duction [ISl, ll9| are applied by reweighting the Z/j* bo- 
son transverse momentum, rapidity, and invariant mass 



distributions from PYTHIA. 

The largest background arises from photon+jets and 
multijet final states in which photons or jets are mis- 
reconstructed as electrons. Smaller background contribu- 
tions arise from electroweak processes that produce two 
real electrons in the final state. The multijet background 
is estimated using collider data by fitting the electron 
isolation distribution in data to the sum of the isolation 
distributions from a pure electron sample and an EM-like 
jet sample. The pure electron sample is obtained by en- 
forcing tighter track matching requirements on the two 
electrons with 80 < Mee < 100 GeV. The EM-like jets 
sample is obtained from a sample where only one good 
EM cluster and one jet are back-to-back in azimuthal 
angle (/). The contamination in the EM-like jets sample 
from W —f ev events is removed by requiring missing 
transverse energy fx < 10 GeV. The average multijet 
background fraction over the entire mass region is found 
to be approximately 0.9%. Other SM backgrounds due 
to W + ^, W+iets, WW, WZ and ti are estimated sep- 
arately for forward and backward events using pythia 
events passed through the GEANT simulation. Higher or- 
der corrections to the pythia leading order (LO) cross 
sections have been applied 19|, l20|, l2]J . These SM back- 
grounds are found to be negligible for almost all mass 
bins. The Z/j* — > t^t~ contribution is similarly negli- 
gible. 

In the SM, the ApB distribution is fully determined 
by the value of sin^ 9^ in a LO prediction for the pro- 
cess qq -> Z/j* -> £+£-. The value of sin2 6'^ is ex- 
tracted from the data by comparing the background- 
subtracted raw ApB distribution with templates corre- 
sponding to different input values of sin^ 9^ generated 
with PYTHIA and GEANT-based MC simulation. Al- 
though sin^ 9^ varies over the full mass range 50 < 
Mee < 500 GeV, it is nearly constant over the range 
70 < Alee < 130 GeV. Over this region, we measure 
sin^ 9^ = 0.2321 ± 0.0018 (stat.) ± 0.0006 (syst.). The 
primary systematic uncertainties are due to the PDFs 
(0.0005) and the EM energy scale and resolution (0.0003). 
We include higher order QCD and electroweak correc- 
tions using the ZGRAD2 [24I program with the generator- 
level Z/j* boson pt distribution tuned to match our 
measured distribution [23| . The effect of higher order cor- 
rections results in a central value of sin^ 9^ = 0.2326 [24|. 

Due to the detector resolution, events may be recon- 
structed in a different mass bin than the one in which 
they were generated. The CC and CE raw ApB distri- 
butions are unfolded separately and then combined. The 
unfolding procedure is based on an iterative application 
of the method of matrix inversion 25| . A response ma- 
trix is computed as Rf^^ for an event that is measured 
as forward in Mee bin i to be found as forward and in 
bin j at the generator level. Likewise, we also calculate 
the response matrices for backward events being found 
as backward {Rf,^), forward as backward {Rf,^), and 



backward as forward {Rf,^). Four matrices are calcu- 
lated from the GEANT MC simulation and used to unfold 
the raw ApB distribution. The method was verified by 
comparing the true and unfolded spectrum generated us- 
ing pseudo-experiments. 

The data are further corrected for acceptance and 
selection efficiency using the GEANT simulation. The 
overall acceptance times efficiency rises from 3.5% for 
50 < Mee < 60 GeV to 21% for 250 < Mee < 500 GeV. 

The electron charge measurement in the central 
calorimeter determines whether an event is forward or 
backward. Any mismeasurement of the charge of the elec- 
tron results in a dilution of ApB- The charge misiden- 
tification rate, /g, is measured using GEANT-simulated 
Z/7* — > e^e~ events tuned to the average rate measured 
in data. The misidentification rate rises from 0.21% at 
50 < Mee < 60 GeV to 0.92% at 250 < Mee < 500 GeV. 
The charge misidentification rate is included as a dilution 
factor V in ApB, with I? = (1 - 2/q)/(1 - 2/q + /|) for 
CC events and V = {1 - 2/q) for CE events. 

The final unfolded ApB distribution using both CC 
and CE events is shown in Fig. [U compared to the 
PYTHIA prediction using the CTEQ6L1 PDFs [ll| and 
the ZGRAD2 prediction using the CTEQ5L PDFs ^. 
The x^/d.o.f. with respect to the pythia prediction is 
16.1/14 for CC, 8.5/14 for CE, and 10.6/14 for CC 
and CE combined. The systematic uncertainties for the 
unfolded ApB distribution arise from the electron en- 
ergy scale and resolution, backgrounds, limited MC sam- 
ples used to calculate the response matrices, acceptance 
and efficiency corrections, charge misidentification and 
PDFs. The unfolded ApB together with the pythia and 
ZGRAD2 predictions for each mass bin can be found in 
Table [D The correlations between invariant mass bins 
are shown in Table HTl 

In conclusion, we have measured the forward-backward 
charge asymmetry for the pp —f Z/j* -I- AT — > e+e~ + X 
process in the dielectron invariant mass range 50 ~ 500 
GeV using 1.1 fb~^ of data collected by the DO exper- 
iment. The measured ApB values are in good agree- 
ment with the SM predictions. We use the ApB mea- 
surements in the range 70 < Alee < 130 GeV to deter- 
mine sin^ 9^ = 0.2326 ± 0.0018 (stat.) ± 0.0006 (syst.). 
The precision of this measurement is comparable to 
that obtained from LEP measurements of the inclusive 
hadronic charge asymmetry [3| and that of NuTeV mea- 
surement y . Our measurements of sin 9^ in a dilepton 
mass region dominated by Z exchange, which is primarily 
sensitive to the vector coupling of the Z to the electron, 
and of ApB over a wider mass region, which is in addition 
sensitive to the couplings of the Z to light quarks, agrees 
well with predictions. With about 8 fb"-'^ of data ex- 
pected by the end of Run II, a combined measurement of 
ApB by the CDF and DO collaborations using electron 
and muon final states could lead to a measurement of 
sin 9'^ with a precision comparable to that of the current 



world average. Further improvements to current MC gen- 
erators, incorporating higher order QCD and electroweak 
corrections, would enable the use of such measurement in 
a global electroweak fit. 
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FIG. 1: Comparison between the unfolded Afb (points) and 
the PYTHIA (solid curve) and ZGRAD2 (dashed line) predic- 
tions. The inner (outer) vertical lines show the statistical 
(total) uncertainty. 



Mee range 
(GeV) 



{Mee) 

(GeV) 



Predicted Afb 

PYTFHA ZGRAD2 



Unfolded Afb 



50- 

60- 

70- 

75- 

81- 

86.5- 

89.5- 

92- 

97- 

105- 

115- 

130- 

180- 

250- 



60 
70 
■75 
81 
86.5 
89.5 
92 
97 
105 
115 
130 
180 
250 
500 



54.5 
64.9 
72.6 
78.3 
84.4 
88.4 
90.9 
93.4 
99.9 
109.1 
121.3 
147.9 
206.4 
310.5 



-0.293 

-0.426 

-0.449 

-0.354 

-0.174 

-0.033 

0.051 

0.127 

0.289 

0.427 

0.526 

0.593 

0.613 

0.616 



-0.307 

-0.431 

-0.452 

-0.354 

-0.166 

-0.031 

0.052 

0.129 

0.296 

0.429 

0.530 

0.603 

0.600 

0.615 



-0.262 

-0.434 

-0.386 

-0.342 

-0.176 

-0.034 

0.048 

0.122 

0.301 

0.416 

0.543 

0.617 

0.594 

0.320 



± 0.066 
± 0.039 
± 0.032 
±0.022 
±0.012 
± 0.007 
± 0.006 
± 0.006 
±0.013 
± 0.030 
± 0.039 
± 0.046 
± 0.085 
±0.150 



± 0.072 
± 0.040 
±0.031 
± 0.022 
±0.014 
± 0.008 
± 0.005 
± 0.007 
±0.015 
± 0.022 
± 0.028 
±0.013 
±0.016 
±0.018 



TABLE I: The first column shows the mass ranges used. The 
second column shows the cross section weighted average of 
the invariant mass in each mass bin derived from pytpha. 
The third and fourth columns show the Afb predictions from 
PYTHIA and ZGRAD2. The last column is the unfolded Afb; 
the first uncertainty is statistical, and the second is system- 
atic. 
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TABLE II: Correlation coefficients between different Mee mass bins. Only half of the symmetric correlation matrix is presented. 



[16] J. Pumplin et ai, JHEP 0207, 012 (2002); D. Stump et 

al, JHEP 0310, 046 (2003). 
[17] R. Brun and F. Carminati, CERN Program Library Long 

Writeup W5013, 1993 (unpubhshed) . 
[18] C. Balazs and C.P. Yuan, Phys. Rev. D 56, 5558 (1997). 
[19] R. Hamberg, W.L. van Neerven, and T. Matsuura, Nucl. 

Phys. B359, 343 (1991); 644, 403(E) (2002). 
[20] J.M. Campbell and R.K. Ellis, Phys. Rev. D 60, 113006 

(1999). 
[21] N. Kidonakis and R. Vogt, Phys. Rev. D 68, 114014 

(2003); M. Cacciari et al, JHEP 04, 68 (2004). 
[22] U. Baur, S. Keller, and W.K. Sakumoto, Phys. Rev. D 



57, 199 (1998); U. Baur et al, Phys. Rev. D 65, 033007 

(2002). 
[23] V. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

100, 102002 (2008). 
[24] This value of sin^ 6^ cannot be compared directly with 

the world average due to the different treatment of elec- 

troweak corrections. 
[25] G.L. Marchuk, Methods of Numerical Mathematics 

(Springer, Berlin, 1975). 
[26] H.L. Lai et ai, Eur. Phys. J. C 12, 375 (2000). 



